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We present a Non-Equilibrium Green's Function (NEGF)-based model for spin torque transfer (STT) devices which 
provides quantitative agreement with experimentally measured (1) differential resistances, (2) Magnetoresistance (MR), 
(3) In-plane torque and (4) out-of-plane torque over a range of bias voltages, using a single set of three adjustable 
parameters. We beheve this is the first theoretical model that is able to cover this diverse range of experiments and a key 
aspect of our model is the inclusion of multiple transverse modes. We also provide a simple explanation for the 
asymmetric bias dependence of the in-plane torque, based on the polarization of the two contacts in energy range of 
transport. 
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Spin Torque Transfer (STT) devices that can switch the 
magnetization of a soft ferromagnetic layer through spin polarized 
electrons without any external field have generated significant 
interest from both basic and applied points of view (See for 
example Ref. [1-15]). Although the concept of spin torque has been 
demonstrated by a number of experiments [14,15], quantitative 
measurement of bias dependence of spin torque has been achieved 
only very recently [16-20]. Note that these experiments show 
considerable dispersions in measuring the bias dependence of in- 
plane torque [16,17,19,20]. Currently, there is no consensus as to a 
microscopic model that accounts for this discrepancy. Moreover, 
the existing theoretical models based on effective mass, tight- 
binding [21-23] and Ab-initio [24,25] band structures do not 
provide quantitative agreements with the experiments. Therefore, 
we need a model which can simultaneously explain all of the 
diverse aspects of STT devices namely i) differential resistances 
R(V), ii) TMR, iii) in-plane/spin-transfer (xy) and iv) out-of- 
plane/field-like (xj components of spin torque. 

This paper presents a simple effective mass model with five 
parameters: a) equilibrium Fermi level Ef, b) Spin-splitting A, c) 
Barrier height of the insulator U/,, d) effective masses for 
electrons inside FM contacts (mpM,f= mpM,j.= mpM) and e) 
effective mass for electrons inside insulator m„^ in terms of which 
we can understand all of the aforementioned characteristics (i-iv) 
of STT devices. Note that we view Ub, mpM and mo^ as parameters 
that account for a wide variety of factors including imperfection at 
ferromagnet/insulator interfaces. As such we consider these three 
parameters adjustable from one structure from another. On the other 
hand, Ef and A are material parameters. Although this is an effective 
mass model that does not include bandstructure effects explicitly, 
we believe that the quantitative agreement with such a diverse set 
of experiments shows that it captures much of the essential physics 
at least in the structures analyzed. One such effect our model tries 
to capture is the role of transverse modes on the TMR, R(V) and 
bias asymmetry of the in-plane torque. This last point is currently a 
topic of debate in the field [16,17,19,20] and our model not only 
fits the data quantitatively, but also provides a natural explanation 
for it. 

Theoretical Formalism: Fig. 1 (a) shows the schematic of a 
magnetic tunnel junction (MTJ) trilayer device that consists of a 
fixed magnet with magnetization M separated by an insulator from 
a free layer with magnetization m, (6 = cos ' (m.M)) being the angle 
between them. Current flows along y, perpendicular to the easy 
plane (x - z) of both the magnets. Our model (See supplementary 
information SI) for spin transport is based on Non-Equilibrium 
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Green's Function (NEGF) [26-28] method using single band 
effective mass Hamiltonian described by the parameters described 
earlier (A, Ub, mpM and m^^ defined in Fig. 1 (b)). The parameters 
Ef, A, and Ub are defined with respect to the bottom of the 
conduction band of the ferromagnet. 




(a) (b) 
Fig. 1(a) Schematic of symmetric tri-layer device. The left contact is the fixed 
FM having magnetization along the z-axis. The right contact is the free layer 
and the channel material is an oxide, z is the easy axis and x-z is the easy plane. 
Transport occurs in y-direction. The device region is modeled using appropriate 
Hamiltonian H and the contacts by self energy matrices Zl,r. (b) The figure 
shows energy-band diagram of the device. The magnets are described by splitted 
bands of up and down spins while oxide only introduces a baiTier. Different 
transverse modes are shown so that conservation of transverse wave vector is 
maintained through the device. 

The choice of the values of Ef and A are independent of the 
specific device because in all the experiments CoFeB is used as 
ferromagnetic contacts. Ef is chosen as 2.25 V, which is believed to 
be typical for CoFeB alloy [19]. The spin splitting A of 3d states of 
transition metals is correlated linearly with the spin angular 
moment mg (See Fig. 3 of [29]). Based on the experimental studies 
involving x-ray absorption and magnetic circular dichroism, mg/usj 
= 0.7 (for Co6oFe2oB2o alloys) [30], where Uj^ the number of d holes. 
Since, exact n^^ of CoFeB is not known, we use the n^^ available for 
bulk Fe (~3.4) to obtain an mg- 2.38 and hence a A ~ 2.38 eV. We 
use a slightly lower value of A = 2.15 eV, which is justified due to 
presence of surface imperfection and other defects at the 
FM/insulator interface. The effective mass mp^ of the CoFeB alloy 
is not reported in literature and is chosen to obtain quantitative 
agreement for the reported differential resistance, R(V) at different 
relative angular positions of ferromagnetic contacts. In practice, the 
values of the effective masses inside ferromagnetic alloys could be 
very different and might depend on the growth of the crystal inside 
the lattice. Among the other parameters, we choose barrier height 
Ub~ 0.7-1 eV and effective mass of electrons inside insulator m^^ ~ 
0.2-0.3 m^, where m^ is the free electron mass. The range of the 
values of Ub is motivated by the presence of oxygen vacancy 
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defects that is expected to reduce the effective Ub of MgO tunnel 
barrier compared to the ideal one (3.7 eV) [31]. The choice of 
also does not differ much from the calculated effective mass of 0.37 
mo from Ab-initio theory [32]. 

An important aspect of our model is the explicit inclusion of 
multiple transverse modes with wave vector k, which are assumed 
to be decoupled like individual 1-D wires. Device Hamiltonian H 
and self energy matrices El,r for each transverse mode are 
described by Eqns. S.I.I and S.I.2. Using H and Zlr, we can 
calculate the Green's Function G(E) and the spin-dependent 
correlation function G'XE), that are used to evaluate charge and spin 
current densities between different lattice sites for each k, (See 
Eqns. S.I. 4, S.I. 5; supplementary information SI). We then sum 
over the results obtained from all transverse modes to get the final 
result relevant to realistic experimental devices having sizable 
cross-sections. 

Note that, in our model we assume that the wave vector kt for 
each transverse mode is conserved throughout the device. The 
transverse energies AEpM.t and AEo^ , in the FM region and the oxide 
(See Fig. 1(b)) are different due to difference in the effective 
masses. In any practical tunnel-junction structure, the barrier is 
unlikely to be defect-free over the junction area because of surface 
imperfections and lattice mismatch between Fe and MgO [32,33]. 
Such effects could lead to non-conservation of kj, but are outside 
the scope of this paper. 
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Experiments 


Ef(eV) 


A(eV) 


mpM /m„ 




m„x /m„ 


Ref. [161, |17] 


2.25 


2.15 


0.73 


0.93 


0.32 


Ref. |19| 


2.25 


2.15 


0.8 


0.77 


0.18 



Table. 1. Parameters Er, A, mpM , Ub and mox used to model two sets of 
available experimental R (V), TMR and Torque. 

R(V) and TMR (Experiment vs. Theory): Most of the 
experiments on MTJs measure differential resistance R(V) that is 
then used to determine TMR (=(Rap-Rp)/Rp). Our model 
successfiilly describes reported R (V), TMR and its roll-off with 
voltage, unlike any other theoretical calculation that we are aware 
of. Fig. 2 (a) and (b) shows comparison between reported [16,17] 
and calculated R (V) at four different relative angular positions (0 = 
0°, 52°, 71°, and 180°) of fixed and free layer magnetizations. The 
parameters are used to benchmark R (V) data with the experiment 
are shown in the top row of Table 1 . 

The simulated bias dependence of R (V) (Fig. 2 (b)) at different 
angles closely agrees with measured data (Fig. 2 (a)). We see that 
for the parallel configuration, Rp has very weak dependence on 
voltage (« 3 kQ). As the angle between the magnetizations of 
ferromagnetic contacts increases, bias dependence of R(V) becomes 
stronger. Inset of Fig. 2 (b) shows bias dependence of TMR. Here, 
the zero voltage TMR (300 K) is calculated to be 150% (indicated 
by 'O')^ which is very close to the reported value of 154%. We also 
calculate the expected roll-off of TMR (43% at 0.54V; indicated by 
'x'). Switching current at parallel configuration is calculated to be 
0.64 mA, which is very close to the reported value i.e. 0.6 ±0.1 
mA. Note that the reported switching current is not directly 
observed, but extrapolated from measured experimental data [16]. 

We also compared simulated differential resistances with the 
reported data of Ref [19] at three different angles G = 0°, 137°, and 
180° (See Fig. 2 (c), (d)). The set of parameters used to benchmark 
R (V) are shown in the bottom row of Table (1). The values of the 
above parameters are chosen based on the arguments described 
above to benchmark the previous experiment. Note that the values 
for all the parameters except Ub and mox remain unchanged. 
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Fig. 2. (a),(c) Experiment and (b),(d) Numerical Simulation of bias dependence 
of differential resistances for different angular positions of free FM w.r.t fixed 
magnet. (Upper panel) [16] and lower panel [19]. Inset of (b) and (d) shows bias 
dependence of TMR, 

Like the previous experiment, bias dependence of R (V) (See Fig. 
2 (d)) closely agrees with those of Fig. 2 (c) for all the angles 
mentioned above. In parallel configuration, similar to the previous 
experiment Rp is roughly constant (~ 122Q). Inset of Fig. 2 (d) 
shows bias dependence of TMR. In this case, we calculate the zero 
voltage TMR (300 K) to be 154% (indicated by 'O'), which is 
exactly same as reported in the experiment. Though the 
experimental switching current density Jc is not reported, based on 
R(V) curve and the switching voltage for P AP transition, we 
can calculate Jc to be 2xl0' A/cm". The numerical value of Jc is 
confirmed by micromagnetic simulation as well as Sun's 
calculation [34] for critical current density. 

Note that we obtain close agreement of R (V) and TMR with the 
measured data from two different experimental groups with a small 
adjustment of Ub and m^^ . 

Comparison of Ty and Tj. (or, 5t||/5V and d^JdW): The quantities 
analyzed above are obtained directly from experiments. However, 
the reported values of Ty and Tx (or, dz\\ldW and dxJdV) are not 
measured directly, rather derived from the measured R.F. voltage 
output (in the spin-transfer-driven ferromagnetic resonance (ST- 
FMR) experiments) using a specific theoretical model with its own 
built-in assumptions. (See Supplementary material [16,17,19]). We 
calculate torques using the same set of values of the parameters 
used to benchmark R (V) with the experiment. 

Fig. 3(a) and (b) show comparison between experimental [15,16] 
and calculated di\\ldW and dxJdW respectively at two different 
relative angular positions (9=58° and 131°) of fixed and free layer 
magnetizations. From our numerical simulation, we see that both 
the torkances are proportional to SinO, so when these quantities are 
normalized with respect to SinO, they should fall onto single curves 
(See Fig. 3 (b)). It confirms the angular dependence of both the 
torkances. Note that the reported data [17] shows 16% variation of 
in-plane torkance for |V|<0.3V) and concluded the earlier measured 
flat nature of (5T||/5V [16] as an artifact of neglecting significant 
terms in the derivation of the RF voltage output in ST-FMR 



experiment. Upper panel of Fig. 3(b) shows the calculated bias 
dependence of in-plane torkance for the bias range of |V|<0.3V, 
which agrees both qualitatively and quantitatively with measured 
data. With regards to the absolute magnitudes of dx\\ldN, our 
simulation shows that near V=0, (5T||/5V)/Sin0 is 0.11(h/2q) kQ'', 
which is in agreement with the reported data i.e. (0.1(h/2q) kQ'') 
[17] and also with the Ab-inito study (0.14(h/2q) kQ"') [23]. In 
contrast to the in-plane component of torkance, the bias dependence 
of the out-of-plane component dxJdV, shows dxJdN'x.V 
(normalized with respect to Sin9), so that = (Aq+AiV^) Sin0. 
Also from the lower panel of Fig. 3(b), it is evident that Ai is 
completely independent of 9 for the tunneling devices. 
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Fig. 3. (a) Measured (Cornell) [17] and (b) simulated bias dependence of in- 
plane and out-of-plane torkances are plotted for two different angles between the 
fixed and soft magnets i.e. 58° and 131°. 
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Fig. 4. (a) Measured (Tsukuba) [19] and (b) simulated bias dependence of 
magnitudes of the simulated in-plane (Spin-Transfer) and out-of-plane (Field- 
like) torques. 

The magnitudes of Xy and x^ reported by Ref [19] is compared 
with the calculated torques as shown in Fig. 4 (a), (b). We use same 
set of values of the parameters used to benchmark calculated R (V) 
data with the experiment. From our simulation, we see that spin 
transfer component of the torque both quantitatively and 



qualitatively matches with the reported value. There is some 
discrepancy in the magnitude of the field like torque Xx, but 
remember that the reported data is based on an indirect 
measurement of torque as reported in literature. Likewise the 
experiment, we plot only the non-equilibrium portion of x^ (i.e. 
AjV). It is important to note that the magnitude of the equilibrium 
part i.e. Aq is comparable to AiV" at lower voltages and can 
manifests itself as an exchange field introducing an asymmetry in 
the R-H loops. 

Now, one can obtain the simulated switching torque amplitudes 
from Fig. 4(a) and (b). Although the switching voltages for both AP 
-> P (-270 mV) and P ^ AP (380 mV) transitions have been 
reported, the reported torque values are available only for AP P 
switching. We have performed micromagnetic simulations to 
confirm that the torque magnitude from the experiment and our 
model do result in switching and are in agreement with the 
threshold values of switching torque [35]. The details of 
micromagnetic simulation are discussed in Supplementary 
information S2. 

Conditions for Asymmetry: As remarked in the introduction, the 
bias asymmetry of Xy (as well as 5X||/5V) has caused much 
discussion in the literature, since many theoretical models [2,24,25] 
predict that X||(V) should be perfectly anti-symmetric (~aV), making 
5X||/5V perfectly anti-symmetric with V, but the results in Fig. (3), 
(4) indicate otherwise. To gain insight, we use our own numerical 
models to look at the quantity (X||(V)+X||(-V))/ (|X||(V)|+|X||(-V)|) for 
different transverse modes individually as shown in Fig. 5(b), 
where the x-axis is the voltage V and the y-axis is A - (Ej-Et). Note 
that (X||(V)+X||(-V))/ (|X||(V)|+|X||(-V)|) is essentially zero as long as 

A-(Ef-E,) >q|V| (1) 

We have found the above criterion to hold for all choices of Ef and 
A so that A-Ef remains constant. One can see from Fig. 5(a) that the 
transverse modes, for which this inequality is satisfied, effectively 
see '"half-metallic" contacts on either side. This is because for any 
transverse mode E,, (Ef-E,) represents the maximum longitudinal 
energy, so that modes satisfying Eqn. (1) (which we call type 'A') 
only one type of spin is injected and detected and it seems that this 
condition leads to perfect anti-symmetry. In an actual device we 
have contributions from both type 'A' modes as well as other 
modes where either one type of spin is injected but both spins are 
detected (type 'B') or both spins are injected and/or, detected (type 
'B" modes). Neither type 'B' nor 'B" mode satisfies Eqn. (1), and 
the degree of asymmetry actually observed in the experimental 
devices suggests that their contribution to Xy dominate that from 
type 'A' mode. 

Reason for Asymmetry: Interestingly, both the experiments and 
our numerical results show that magnitude of the torque X||(V) is 
smaller when a negative voltage is applied to the free layer (or, 
positive voltage to fixed layer as shown upper panel of Fig. 4) so 
that free layer bands are raised up. Asymmetry in the magnitude of 
X|| is explained by Fig. 5(a). Let us assume the fixed layer is along 
X -direction, while free layer is along z -direction. Therefore, 
according to Eqn. (S.I.6) {supplementary information SI), in-plane 
torque on the free layer is X||_ a-ee (per Bohr Magneton) = l^z, while 
torque acting on the fixed layer is Xy fj^eii (per Bohr Magneton) = - 
Isx- Now, applying a negative voltage to free layer is the same as 
applying positive voltage to fixed layer and we expect torque 
exerted on the free layer at negative voltage to be the same as that 
exerted on the fixed layer at positive voltage. This is indeed true as 
evident from the Isz and Isx plots shown in Fig. 5(c). 
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Fig. 5. (a) The schematic shows the fixed and free layers are along z and x-axis 
respectively. Ferromagnet band-structure is drawn explaining the relative 
contribution of transverse modes on the asymmetric bias dependence of i||. (E|- 
E,) represents maximum longitudinal energy for which there are thi'ee types of 
modes: type (A) (E|-E,<A-qV), type (B) (A-qV< E|-Ei<A), and (B') modes (E,- 
E,>A). ErEi vs. Pc is drawn for both the FM contacts showing relative 
contributions of density of states of majority and minority earners available for 
tunnehng. Parameters used to generate (b), (c) and (d) are same as shown in the 
top row of Table 1. (b) shows the plot of the fUnction (x||(V)+X||(-V))/ 
(|t||(V)|+|t||(-V)|) for different transverse modes, the x and y axes are being 
voltage V, and A-(E|-E,). Black denotes a value of the fiinction less than a part 
per milhon (10"'). It separates the regions of anti-symmetric and asymmetric 
bias dependence of ly. The two variables V and A- (E| -E,) have been varied 
independently to produce 2D grid with 2,50,000 grid points, (c) Spin cuirent is 
calculated for type (B') modes; When negative voltage is applied to fixed layer, 
we get Is£>Isx. When the voltage polarity is reversed, we get Is£<Isx. Bias 
dependence of spin currents for type (B) modes is same as that of (B') modes, 
(d) Spin current is calculated for type (A) modes. We see that \sz " Isx 
irrespective of voltage polarity. We see that the magnitude of spin current is 
much greater for both type (B) and (B') modes compared to type (A) modes. 
Therefore, type (B) and (B') modes contribute more in tunneling compared to 
type (A) modes for the devices explored in the literatures. 

Note that Ig^ > Ig^ (for type 'B" (as well as 'B') modes), when 
negative voltage is applied to fixed (z) magnet w.r.t. free (x) 
magnet (Fig. 5(c)). But Isz = hx for type 'A' modes. This can be 
understood by considering spin polarization (Pc) of the fixed (z) 
magnet and free (x) magnets in the energy range where current 
flow occurs for a given transverse mode (See Fig. 5(a)). 
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Energy range 
of transport 


Pc(i) 


Pc(X) 


^11 (V) 




Ef-E, < A-qV 


1 
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Anti-symmetric 
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A-qV <Ef-E, <A 


1 


<1 


Asymmetric 


B' 


Ef-E,>A-qV 


<1 


<1 


Asymmetric 



Table. 2. The modes coiTesponding to particular energy range of transport and 
bias dependence of in-plane torque,!]. 

When the polarity of the voltage is reversed, for the type 'B" 
modes, Pc (z) <Pc (x) <1 and for the 'B' modes, Pc (z) <Pc (x) = 
1. In both cases, Isz < Isx (See Fig. 5(c)). For type 'A' modes, Pc (z) 



= Pc (x) =1 and hence, Isz=Isx (Fig. 5(d)). Note that the out-of- 
plane component is determined by Isy whose magnitude depends on 
the Pc of both the contacts and hence unaffected by the voltage 
polarity. By contrast the in-plane torque on the free layer is 
determined by Isz whose magnitude depends on the Pc of the fixed 
(z) magnet and is larger when a negative voltage is applied to it. 
This seems to provide a simple explanation for the asymmetry of 
the in-plane torque (top panel of Fig. 3, 4) which attracted much 
discussion. 

To summarize, we have described a quantum transport model for 
a spin torque device using NEGF formalism, which provides both 
qualitative and quantitative agreement on differential resistances 
R(V), TMR, and torkances simultaneously with a single set of 
parameters (i.e. effective masses, barrier height and spin-splitting in 
contacts). We show that bias asymmetry in Xy can be understood in 
terms of the polarization (Pc) of the fixed layer in the energy range 
of transport. Though our analysis is focused on vertical STT 
devices, the general formalism can also be used for analyzing bias 
dependence of spin torque in non-local lateral spin-torque devices 
[36,37]. 

The work is supported by MARCO focus center for Materials, 
Structure, and Devices. 
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Supplementary Information 

Voltage Asymmetry of Spin-Transfer Torques 

Deepanjan Datta^, Behtash Behin-Aein^, Sayeef Salahuddin^'^, and Supriyo Datta'^ 
School of Electrical & Computer Engineering, Purdue University, West Lafayette, IN-47907, USA ^ 
Department of Electrical Engineering and Computer Science, University of California, Berkeley, CA-94720, USA 

SI. Details of Non-Equilibrium Green's Function (NEGF) Model used 

The Non Equilibrium Green's Function (NEGF) based transport model for MX J devices is based on 
single band effective mass Hamiltonian H and self-energy S^r which are used to calculate Green's 
fimction G(E), electron correlation fimction G°, charge and spin current densities J and f' respectively. 

(a) Hamiltonian and Self-Energy matrix 

Here, we will present the device Hamiltonian H and self energy matrices Sl,r for each transverse 
mode kj through the device. We assume that the fixed magnet M and the soft magnet m are both in the 
X - z plane (see Figure SI). Hamiltonian for each transverse mode is given by (ct = Pauli Spin matrices) 



Left Contact: Hl(?,7, kj = 



(a,^(k.)+ qV/2) I 



+ 



^I-a.M^ 



A, i = j 



- tp]^ I , if i and j are nearest neighbors 



0, otherwise 



(S.Ll.a) 



Left Interface: H,„terface(^^^^k.) = ("..(k.) + (Ub/2) + (qV/2) )l + 



I - a.M 



A (S.I.I, b) 



Channel: 



Hcliamel(^'./' kt) " 



aox(kt) + U,+qV(i--^))l, i = j 
X^^ I , if and j are nearest neighbors 



0, otherwise 



(S.I.1.C) 



Right Interface: = (aM+ (Ub/2) - (qV/2))l + 



I - a.m 



(S.I.I, d) 
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Right Contact: H,^ k, ) = 



(a,^(k,)-qV/2)l 



+ 



- a.m ^ 



A, i = j 



I , if i and j are nearest neighbors 



0, otherwise 



(S.I.l.e) 



where, (k, ) = 2t^ + hX I , (k, ) = 2t„, + /2m*^ . {K ) = 0.5 (a^^, + ) , 
tpM ~ ^^/2nip]y[a^ , t^^ = h^/2m^^a^ , nipj^^and m*^are effective masses of electron inside FM and 

insulator region respectively and a is the uniform lattice spacing and I is the 2x2 identity matrix. N is 
the total number of atomic sites (excluding the interface points) inside the insulator. V is the applied 
voltage difference between two FM contacts which is assumed to drop linearly inside the insulator. is 
the barrier height of the insulator. 

To write the self-energy matrices Zl,r we note that their non-zero elements are the left (L) or, right 
(R) lattice points respectively and these are given by a 2x2 matrix of the form 



FMexp{ik^^a) 








■tFMexp(iki;^a) 



(S.I.2. a) 



where. 



kLRa = cos"^ 



^ E + qV/2 - (h^k?/2m;^) 
2t 



(S.I.2. b) 



1 i -1 
k^^a = cos 



E+ qV/2-(hX/2m;^)-A ' 
2t, 



"-FM 



(S.I.2. c) 



(Use +qV/2 for the contact L, and - qV/2 for the contact R) and E is energy. The form given above for 

2^L,R ) is appropriate, if the easy axis of the magnet is used as the quantization axis for the spin. If 

the easy axis does not lie along z , then a unitary transformation is needed to transform this quantization 
axis to z [SI]. 



6 




Figure SI. (a) The schematic of FM/Insulator/FM device is shown. The device region is modeled using 
appropriate Hamiltonian [H] and the contacts are taken into account by self energy matrices Sl,r whose 
anti-Hermitian components rL,R = i(SL,R -2l,r') are broadening matrices due to contacts L and R 
respectively. ^'^ (b) tpM and tox represent the coupling parameter between each lattice site given by tpM = 
h^/2mFM*a^, tox = hV2mox*a^, a being uniform lattice spacing in discrete representation. apM and aox are 
given by apM = 2 tpM + h\,V2mFM* and aox = 2 tox + h\t^/2mox*, where, h\t^/2mFM* and h\^/2mox* 
represent transverse energy of a mode with wave vector kt inside FM and insulator regions respectively. 
Jspin Left and J Spin Right arc current densities going into and coming out of the soft ferromagnetic region. 
Therefore, the current density absorbed is Jspin fm = Jspin Lefi - Jspin Rigiit, which exerts torque on the soft FM 
layer. 
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(b) Green's Function, Charge and Spin Current Density 

In Appendix Sl(a), we showed the device Hamiltonian H and self-energy matrices SL,Rfor the device. 
Once H and IiL,R are known, all quantities of interest including Green's function G(E) can be calculated 
from the following set of equations [S2, S3] 

Green's Function: G(E) = [E - H - Sl - Sr] ; 

Spectral Density: A = i (G - ) = G T G^; 

Electron Correlation Function: G" (E) = G(T}1 + Sr ) G^; j 3) 

In-scattering Function: r (E) =rL j;^(E) fL j;^(E); 

Broadening Matrix: rLj^(E) = I^Zl^r (E)-Zl,r (E)); 

where, G° {= -iG^) refers to spin-dependent correlation function whose diagonal elements are electron 
density. A is the spectral function, whose diagonal elements are local density of states, is the in- 
scattering function representing the rate at which electrons come in to the device from the contacts, and 
fL,R are the Fermi functions in contacts 1 and 2 respectively. Once we know G", the charge and spin 
current density between lattice sites are calculated as [S3] 



Charge Current Density: j+i = 7^ j dE Real Trace (H^. ^.^jG"+j ^. - G"^.+jH].^j ^. ] 



Spin Current Density: fj j^^ ^ ^ j Real Trace {d. (h^. ^.^iG°^i ^. - G°^^.^iHJ^i ^. )} 

1^ TJ ~ 



(S.I.4) 



(S.I.5) 



The spin torque exerted on the free magnetic layer will be the difference between spin current going into 
and coming out of the soft magnet in the absence of scattering. We define spin torque as the divergence of 
spin current [S4] 

T = - ^3 jdF [v.Js] = la JdS (-jdy V.Js) = ^3 jdS Jspin.Right) = ^^efdS (S.I.6) 

£i S S S 

where, S and V are the area of cross-section and the volume of magnet respectively, i^b is the Bohr 
Magneton (i.e. magnetic moment per unit volume) and Js is the spin current density between two lattice 
points. In our calculations, however, |dS Jgpin fm is replaced by |dS Jgpij, Leg based on following two 
assumptions: (i) Spin-polarized electrons becomes completely polarized along the free layer 
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magnetization m after traversing a few monolayer inside the magnet, (ii) Spin torque acting on the free 
layer has components only transverse to the free layer magnetization. 
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S2. Micromagnetic Simulation 



Current induced magnetization switching in Spin Torque Transfer (STT) devices is of great interests 
because of its range of potential applications. To integrate the STT devices with existing semiconductor 
technology, lower switching current densities are very important. The current density for switching the 
magnetization of the free layer is derived by Sun [S5] 



a 



(tpMHcMs) 



1 + 



27iM, H 

- + - 



H 



The corresponding condition for threshold value for the switching torque is [S5] 



H 



S-min 



= a 



1+ 



- + - 



H 



(S.II.1) 



(S.11.2) 



where, Hs,min is the minimum field required for switching (due to spin torque), He is the uniaxial 
anisotropy field, and a is the damping parameter. The effect of demagnetizing field (Hd) is included by 
the term IttMs/Hc. H includes the external magnetic field and the field due to out-of-plane component of 
spin torque. tpM is the thickness of the thickness of the free FM layer. Now, and demagnetizing field 
of the free layer are 45 Oe and 1.38 T respectively [86]. Therefore, with applying any external field and 
assuming the damping coefficient a to be 0.008, the condition for svwtching becomes Hs,mm = 55.8 - 

0.008 Hout-of-plane- 

Now, we will describe the magnetization dynamics with macrospin approximation. If spin torque T is 
exerted on the free layer vwth magnetization in, the rate of change of m is described by landau-Lifshitz- 
Gilbert's Eqn. and is written as. 



dm I I " , dm ^ 

=- Y mxH fl. +amx +x 

dt ' ' dt 



(S.II.3) 



where, H_.^ (^=H^^j+H(,+H^^ can be vratten as the gradient of the potential energy w.r.t. normalized 
magnetization i.e. Yi^ = -(l/MgF) V^^E where Ms is the saturation magnetization. The above equation 
can be simphfied as. 



dm 
"dT 



H 



a ^ 



eff 



(S.II.4) 



We include the effect of spin torque by defining the components of xas 
T^=T^(V)(mxM^ and\=x^{y) mx^mxM^ . Substituting in Eqn. (S.II.4), we get 
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(l +a')-^ =-|y| mx(H^g+AM ) -a|Y| mx^mx^H^^+B M )) (S.II.5) 

where, the following notations have been adopted: 

A= T^(V) -aT^(V) B= T^(V) + T^(V)/a 

Now, we solved the Eqn. (S.II.5) with varying the magnitude of X||(V) and t_l (V) calculated for different 
voltages. We see that AP ^ P switching occurs when a negative voltage is apphed to the fixed layer (V = 
-0.27 V) which is exactly same as the switching voltage mentioned in Ref. [S6]. We see that at V = -0.27 
V, magnitude of the in-plane torque is H||(V=-0.27V) = T||(V=-0.27V)/MsV = 60 Oe and H^(V=-0.27V) = 
Tx (V=-0.27V)/MsV = 40 Oe. Now if we put the values of in Eqn. (S.III.2), Hs,„^= 55.8 Oe < H||(V=- 
0.27V). Therefore, the condition for switching derived by Sun is satisfied with the values we obtained 
from NEGF simulation. Figure S2 shows the magnetization vs. time diagram of our micromagnetic 
simulation. 
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Time (nS) 



Figure S2. Micromagnetic simulation for AP P switching. The parameters used are: He = 45 Oe, 
H||(V=-0.27V) = T||(V=-0.27V)/MsV = 60 Oe and Hi(V=-0.27V) = (V=-0.27V)/MsV = 40 Oe, Mg = 
1 100 emu/ cm^, tpM = 2 nm and a = 0.008. The switching occurs exactly at V = -0.27 V which is reported 
in paper. For this simulation we have taken into account the out-of-plane demagnetizing field Hd of 
1.38 T. 
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